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ABSTRACT 

We derive the Y-D relation of Galactic supernova remnants of shell-type separately at 
adiabatic-phase and at radiative-phase through two sets of different formulas, considering the 
different physical processes of shell-type remnants at both stages. Also statistics on Galactic 
shell-type remnants about 57 was made. Then we do some comparison with other results ob- 
tained before. It shows that all the best fit lines in the Y-D relation plots newly are to some 
extent flatter than those derived by some authors at early time. Our theoretical and statistical 
outcomes are in somewhat good consistency. 

Subject headings: methods: — statistical — analysis — (ISM:) supernova remnants 



1. Introduction 

The relation between the radio surface bright- 
ness (E) and the diameter (D) of SNRs has being 
widely discussed before (e.g., Poveda & Woltjer 
1968; Clark & Caswell 1976; Green 1984; Mills et 
al. 1984; Huang & Thaddeus 1985; Arbutina et al. 
2004, etc.). Many authors use them to determine 
the distance of a SNR (Poveda & Woltjer 1968; 
Clark & Caswell 1976; Lozinskaya 1981; Huang & 
Thaddeus 1985; Duric & Seaquist 1986; Guseinov 
et al. 2003). There are some different outcomes 
about the best fit line of the E-D relation. De- 
spite that one straight line was derived by some 
workers in their statistics (e.g., Poveda & Woltjer 
1968; Huang & Thaddeus 1985; Arbutina et al. 
2004), Clark & Caswell (1976), Allakhverdiyev et 
al. (1983), and Allakhverdiyev et al. (1985) got 
a broken fit line in their previous derivation. At 
408 MHz Clark & Caswell (1976) have a broken 
line with slopes of /3 = -2.7/ — 10 (E oc Z^)at 
D < 32 pc/D > 32 pc. While at 1 GHz Al- 
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lakhverdiyev et al. (1985) gave line slope of —3.0 
and —6.0 with broken point at D = 40 pc. The- 
oretical analysis once made by Duric & Seaquist 
(1986) got following results 



Y(D) = 4 x lO" 15 !)- 3 - 5 , D < lpc 
Y(D) = 4 x 10- 14 L>- 5 , D > lpc 



(1) 
(2) 



With the line broken at 1 pc when the remnants 
are far too small to be detected. 

Galactic SNRs are classified into three types: 
Shell-type, Plerion-type and Composite-type. 
Merely shell-type remnants are analyzed in our 
work. Furthermore, shell-type SNRs usually 
have four evolution stages: the free expansion 
phase, the Sedov or adiabatic phase, the radia- 
tive or snowplough phase and the dissipation 
phase. Nearly all of the detected SNRs are in 
the adiabatic-phase, or in the 3rd. Almost none 
is observed in the 1st and 4th phases. Here we do 
a simple statistics on the Y-D relations of some 
Galactic SNRs in section [2j and do some theoret- 
ical reductions of Y-D relations independently at 
adiabatic-phase and radiative-phase in section [3J 
and some comments made in section [4] In the last 
section summarizes our conclusion. 
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2. Statistics of E-D relation 

We newly collected 57 shell- type remnants (ta- 
ble[5]) in Galaxy with known distance (d) which are 
derived not by E-D relation but by some other dif- 
ferent methods to make statistics, and get relation 
as follow, 

53(D) = 5.50 x 10- 18 D- 215±o - 38 

(Wm' 2 Hz~ 1 sr~ 1 ) (3) 

The SNR G4.5+6.8 (Kepler) in table [5] was ex- 
cluded since its extremely small linear diameter 
(3 pc only) so as to avoid large deviation of fit- 
ting. Another source SNR G166.2+2.5 was also 
not used for being a false remnant. In Fig. [I] the 
E-D relation was shown. It seems that the best 
fit lines in plot are somewhat flatter than those 
derived by some authors at early time (e.g., Clark 
& Caswell 1976, 1979; Milne 1979; Allakhverdiyev 
et al. 1985; Duric & Seaquist 1986). In the follow- 
ing we do some theoretical works about the E-D 
relation of the remnants. 
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Fig. 1. — For 57 Galactic SNRs of S-type the sur- 
face brightness (E) decreases with their linear di- 
ameter (D) by a slope of about —2.15 ± 0.38. 

3. Theoretical reduction 

The physical processes of supernova remnants 
at Sedov-phase are rather different with that 
at radiative-phase. Therefore we do theoretical 
derivation of the E-D relation of remnants sep- 
arately for both stages, since SNR at both free- 



expansion phase and last phase are practically 
undetectable. 

3.1. E-D relation at Sedov-phase 

Assuming the remnants linear diameter (D) in 
pc, SNe initial explosion energy (D 51 ) in the unit 
of 10 51 ergs, and ISM electron density (n ) in 
cm~ 3 , from the standard Sedov solution we have 
the following well-known equation for remnants at 
second stage (Itoh 1978, Bignami & Caraveo 1988, 
Zaninetti 2000, Volk et al. 2002, Ptuskin & Zi- 
rakashvili 2003) 

D(t) = 4.3 x HT 11 ( ) V5 i 2/5 (4) 

After making differentiation dD Jf^ (Duric & 
Seaquist 1986), one has 

v{t) = lA Q t-^ 5 (5) 
5 

where 

A = 4.3 x HT 11 ( f ) 1/5 (6) 

At adiabatic phase the remnants thicknesses are 
proportional to D (Milne 1970). Then we have 
the shell volume as 

V(D) = C D 3 (7) 

Here 

c » = iOU) 3 H- 37 < 8 » 

when we approximately assume Di/D Q ~ 2/3. 
Di(D ) is the inner(outer) diameter of the rem- 
nant shell. Combining (3) and (6), one gets 

V(t) = CoAlt 6 / 5 (9) 

As the remnants shock waves travel at the Sedov- 
phase, the ambient magnetic field B will decrease 
with D according to (Padmanabhan 2000) 

B(D)=B Q {^) 2 (10) 
Substituting (3) to it, we have 

B{t) = B D 2 A 2 t- 4 / 5 (11) 
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While Duric (2000) shows that at the adiabatic 
phase of remnants evolution, their radio emissivity 
e{B,v) is expressed as 



e(B,v) = 3 x 10~ 34 n ( 



B 



7Q00Kms- 



10- 4 G / 
\ GHz J 

{WHz^rn 



(12) 



Where, the magnetic field B is expressed in units 
of 10~ 4 G, the velocity of shock wave in units of 
7000 km s _1 , and radiation frequency in GHz. by 
(4), (10) and taking the average value of the rem- 
nants spectral index a — 0.5, we can get 

e(t) = 0.1714 xHT 4 VD ( f' 2 (13) 

when the shell volume is taken to be the radiat- 
ing electrons encompassing volume, the remnants 
surface brightness can be given by 



E(t) = 



0.1714 x 10 



<t)V{t) 
n 2 D 2 (t) 

-40 n oA) 



x ( T?F*G ) CW~ 3/5 



(14) 



combining (3) and above formula, one gets the fi- 
nal form 



£(£>) = m a D~ 15 (Wm~ 2 Hz~ x sr^) (15) 



where 



0.1714 x 10 



-40 



W- A G 



CoAo - 5 



7.56 x 10~ 19 



(16) 



Here the typical value of some SNRs physical pa- 
rameters is taken: ISM density no = 0.1 cm~ 3 , 
SNe initial explosion energy E§i = 1, the remnants 
diameter and ISM magnetic field at the beginning 
of Sedov phase D — 2 pc and B = 10~ 4 G, etc. 

The derived line in Y.-D relation plot of the 
SNRs in second phase is somewhat flat 



E(D) = 7.56 x lO" 19 !}- 1 ' 5 



(17) 



3.2. Yi-D relation at radiative-phase 

Similarly assuming the remnants linear diame- 
ter (D) in pc, SNe initial explosion energy (£51) 
in the unit of 10 51 ergs, and ISM electron density 
(no) in cm -3 , we have the following equations for 
remnants at third stage (Itoh 1978) 



D(t) = 1.42 ( ) 



5/21 



2/7 



(18) 



Which is rather deferent from formula (3). After 
making differentiation dD J^ , we have 



dt ' 



,{t) = L -Aot-^n 



where 



^0 = 1.42 Sll 



5/21 



(19) 
(20) 



the same as that in Sedov phase 

V(D) = C D 3 (21) 
when we approximately take Di/D D ~ 3/4, then 



Co = - 6 [ 1- 



* 1 < D t " 3 



0.3 



(22) 



Here, Di(D Q ) is defined as before. Changing the 
variant D to t, one have the form 



V(t) = CoAlf' 7 



(23) 



Because that as the remnants shock waves 
travel at Sedov-phase the ambient magnetic field 
B of a remnant decreases with the diameter D ac- 
cording to ([TO]) and at dissipation-phase according 
to B(D) = B (Do/D)°, one can moderately sup- 
pose that at radiative-phase the magnetic field B 
will decrease with D as follow 



B(D) = ^B 



(24) 



which obviously is not the same as the formula (9). 
After substituting (7) to it, one has 



B ^ = ( ^? ) ^ 2/? 



(25) 



At the radiative phase of remnants evolution we 
could assume B ~ 10~ 6 Gauss, and vo ~ 220 Km 
s _1 for the formula in Duric (2000) 



n a+1 



e(t) = 3 x 10~ 34 noA) ( y&a ) 

xA - a - H -H^) ( p 

x ( mn P (26) 
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Where, the magnetic field B is expressed in units 
of 1CP 6 G, the velocity of shock wave in units 
of 220 km s _1 , and radiation frequency in GHz. 
Therefore the remnants surface brightness is ex- 
pressed as 



£(t) 



<t)V(t) 
n 2 D 2 (t) 

3xl0 -3 4( r } « ( )D . + i 



220/f rns" 1 / 



xA^ 3a Co( 

x^P^ 77 (27) 
changing the variant t to D, we have the form 

£(£>) =m r D- 2 (Wm^Hz^sr' 1 ) (28) 
where 

m. r = 3xl0- 34 ( 1 ) 2 «5( ^ 



xC Ao ba = 4.63 x l(T i8 (29) 

Here we have used the typical values for the phys- 
ical initial parameters as that in section f3. 21 

Thus a little steeper line of £-£> relation of 
SNRs at 3rd phase is obtained. 



E(D) = 4.63 x lCT^-Cr 2 



(30) 



3.3. Transition from Sedov-phase to radiative- 
phase 

Assuming D t the diameter of remnant when 
SNR physical state transfers from Sedov phase to 
radiative phase, according to formulae (14) and 
(27), we know 



(31) 



Thus, one gets D t ~ 38 pc which is also shown in 

Fig. ® 

Through statistics, Clark & Caswell (1976) got 
the diameter at broken point about 32 pc for 
29 galactic SNRs at 408 MHz, and also 32 pc 
at 5000 MHz. Allakhverdiyev et al. (1983) got 
D t ~ 30 pc at 408 MHz for 15 shell-type remnants, 
and ~ 32 pc at 1 GHz. For a larger number of 
samples of 146 all-sort galactic objects including 
plerion, shell and composite-type remnants, Al- 
lakhverdiyev et al. (1985) obtained D t ~ 40 pc at 



1 GHz. Our result by theoretical analysis method 
is in somewhat good consistency with theirs. But 
some authors plotted only a straight line by their 
statistics early (Duric & Seaquist 1986; Guseinov 
et al. 2003). 

Obviously the D t value is the average one of 
the diameter of all galactic SNRs. When the den- 
sity (no) of interstellar media or magnetic field 
strength (B) is somewhat larger than usual, or 
so, then the transformation diameter D t may be 
smaller, and vice-versa. 
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Fig. 2. — Comparison plots of some results de- 
rived from Case & Bhattacharya (1998), Xu et 
al. (2005) and on theoretical and statistical 
ones in this paper. The best fit values are 
f3 = -2.64, -1.6, -1.5/ - 2.0 (for Sedov/radiative 
phase) and —2.15 (E oc D@), respectively. All of 
the fit lines are somewhat flatter than those de- 
rived by some authors at early time. 



4. Discussion 

4.1. Results comparison 

Comparison of some results derived from Case 
k Bhattacharya (1998), Xu et al. (2005) and 
on theoretical and statistical ones in this paper 
was made in Fig. [2] The best fit values are 
f3 = -2.64, -1.6, -1.5/ - 2.0 (for Sedov/radiative 
phase) and —2.15 (E oc D 13 ), respectively. All of 
the fit lines are somewhat flatter than those de- 
rived by some authors at early time. 

Case & Bhattacharya (1998) got one straight 
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slightly steeper line 



S(D) = 5.43 x 1Q- 17 D- 2 M 



pc 

(Wm^Hz^sr- 1 ) (32) 
And Xu et al. (2005) derived a straight flatter line 



E(L>) = 1.21 x 10- 18 £>- 160 



(Wm- 2 Hz- 1 sr- 1 ) 



(33) 



The best fit line also shows no broken in our new 
statistics of this paper. Broken in our theoretical 
analysis is to some extent very small. 

4.2. SNR total fluxes increase with age at 
Sedov-phase? 

The radio surface brightness of a remnant de- 
fined as 

s = Sum, (34) 

Here, S\ghz is the detected flux of a remnant at 
1 GHz, 9 is the observational angle. 
And we know 



9 2 X D 2 

Therefore one has 

E oc SiGHzD~ 2 

Let us suppose that 

SlGHz = So,igh z D^ 



(35) 
(36) 
(37) 



Where, Sq,ighz is the SNR original flux at 1 GHz, 
(3' is an index of the formula above. 
Then 

E cx LH 2 -/ 3 ') = D-P (38) 
Namely (3 = 2-/3' Now we can see that 

1. If f3' = 0, then E cx D~ 2 , (3 = 2, S 1GHz = 
Sa.iGHz and E oc D~ 2 . It means that the 
radio flux of a evolved SNR remains a con- 
stance. 

2. If (3' < 0, then (3 > 2, S 1GHz < S ,igh z - 
In the case, the flux of a evolved remnant 
decreases with age. 

3. But if (3' > 0, then (3 < 2, S 1GHz > S ,igh z - 
the evolved remnant flux increases against 
time. 



Our theoretical outcome and a statistical result 
(Xu ct al. 2005) (/?' = 0.50 and 0.40) are belong to 
the last case, i.e. the SNR radio flux increases with 
age. But our statistics in the paper shows (3' < 0. 
In summary, we are not very certain whether the 
SNR total fluxes increase or decrease with age at 
Sedov-phase. In other words, the S-D relation is 
not significantly sensitive to be used to determine 
the variety of SNR radio fluxes. 

5. Conclussion 

From the basic publicly accepted formulae and 
some simple reasonable physical suggestions, we 
derived mathematically the relation between sur- 
face brightness and linear diameter of Galactic 
shell-type SNRs. Firstly, we did a simple statistics 
about the E-Z) relation of 57 Galactic remnants. 
The diameter at phase-transition from adiabatic- 
stage to radiative-stage is also arithmetically ob- 
tained and ultimately in agreement with the sta- 
tistical results gained before. The line slope val- 
ues in theoretical relation plot are (3 = —1.5/ — 2.0 
for Sedov/radiative phase (E cx D 13 ) respectively. 
Different statistical results derived now and before 
were compared with our theoretical one. It shows 
that all the best fit lines in plots are somewhat 
flatter than those derived by some authors earlier. 
One can guess that even the undertook samples of 
galactic SNRs increase in the future, our statisti- 
cal outcome of the E-D relation will alter a rather 
little. And also our theoretical result on the rela- 
tion can by and large remain valuable. Moreover, 
newly better methods to determine the distance of 
a common supernova remnant are greatly needed. 

JWX likes to thank J. S. Deng and Y. Z. Ma for 
their assistance and help during the paper work. 
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Table 1 

Some physical parameters of 57 shell-type Galactic SNRs. 



Source 


Age 


JJ1SL. 


uia. 


size 


c 

JIGHz 


xtei. 




yr 


pc 


pc 


arcniin 


Jy 




G4.5+6.8 


380 


2900 


3 


3 


19 


H90, G04a 


G7.7-3.7 


— 


4500 


29 


22 


11 


M86 


G8. 7-0.1 


15800 


3900 


51 


45 


80 


G96 


G18.8+0.3 


16000 


14000 


57 


17x11 


33 


D99,G04a 


G27.4+0.0 


2700 


6800 


8 


4 


6 


C82, G04a 


G31.9+0.0 


4500 


7200 


13 


7x5 


24 


CS01 


G32.8-0.1 


— 


7100 


35 


17 


11 


K98b 


G33.6+0.1 


9000 


7800 


23 


10 


22 


S03, SV95, G04a 


G39. 2-0.3 


1000 


11000 


22 


8x6 


18 


C82 


G41. 1-0.3 


1400 


8000 


8 


4.5x2.5 


22 


C82, B82, C99 


G43.3-0.2 


3000 


10000 


10 


4x3 


38 


L01 


G49.2-0.7 


30000 


6000 


52 


30 


160 


KKS95, G04a 


G53.6-2.2 


15000 


2800 


24 


33x28 


8 


S95, G04a 


G55.0+0.3 


1100000 


14000 


71 


20x15? 


0.5 


MWT98 


G65.3+5.7 


14000 


1000 


78 


310x240 


52 


LRH80, R81 


G73.9+0.9 


10000 


1300 


8 


22? 


9 


L89, LLC98 


G74.0-8.5 


14000 


400 


23 


230x160 


210 


LGS99, SI01, G04a 


G78. 2+2.1 


50000 


1500 


26 


60 


340 


LLC98, KH91 


G84.2-0.8 


11000 


4500 


23 


20x16 


11 


MS80, M77, G04a 


G89.0+4.7 


19000 


800 


24 


120x90 


220 


LA96 


G93.3+6.9 


5000 


2200 


15 


27x20 


9 


L99, G04a 


G93.7-0.2 


— 


1500 


35 


80 


65 


UKB02 


G109. 1-1.0 


17000 


3000 


24 


28 


20 


FH95, HHv81, G04a 


Gill. 7-2.1 


320 


3400 


5 


5 


2720 


TFvOl 


G114.3+0.3 


41000 


700 


15 


90x55 


6 


MBP02, G04a 


G116. 5+1.1 


280000 


1600 


32 


80x60 


11 


RB81, G04a 


G116.9+0.2 


44000 


1600 


16 


34 


9 


KH91, G04a 


G119.5+10.2 


24500 


1400 


37 


90? 


36 


MOO 


G120. 1+1.4 


410 


2300 


5 


8 


56 


H90, G04a 


G127.1+0.5 


85000 


5250 


69 


45 


13 


FRS84 


G132. 7+1.3 


21000 


2200 


51 


80 


45 


GTG80, G04a 


G156.2+5.7 


26000 


2000 


64 


110 


5 


RFA92 


G160.9+2.6 


7700 


1000 


38 


140x120 


110 


LA95 


G166.0+4.3 


81000 


4500 


57 


55x35 


7 


L89, KH91, G04a 


G166.2+2.5 


150000 


8000 


186 


90x70 


11 


RLV86 


G182.4+4.3 


3800 


3000 


44 


50 


1.2 


KFR98 


G205.5+0.5 


50000 


1600 


102 


220 


160 


CB99 


G206.9+2.3 


60000 


7000 


102 


60x40 


6 


L86 


G260.4-3.4 


3400 


2200 


35 


60x50 


130 


B94, RG81 


G266.2-1.2 


680 


1500 


52 


120 


50 


K02, AIS99 


G272.2-3.2 


6000 


1800 


8 


15? 


0.4 


D97 


G284.3-1.8 


10000 


2900 


20 


24? 


11 


RM86 
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Table 1 — Continued 



Source 




Age a 


Dist. a 


Dia. b 


size 


SlGHz 


Ref. 






yr 


pc 


pc 


arcmin 


Jy 




G296.5- 


flO.O 


20000 


2000 


44 


90x65 


48 


MLT88 


G296.8- 


-0.3 


1600000 


9600 


47 


20x14 


9 


GJ95, G04a 


G299.2- 


-2.9 


5000 


500 


2 


18x11 


0.5 


SVH96 


G309.2- 


-0.6 


2500 


4000 


16 


15x12 


7 


RHS01 


G315.4- 


-2.3 


2000 


2300 


28 


42 


49 


DSM01, G04a 


G321.9- 


-0.3 


200000 


9000 


70 


28 


13 


SFS89, S89 


G327.4- 


f0.4 




4800 


29 


21 


30 


SKR96, WS88, G04a 


G327.6- 


fl4.6 


980 


2200 


19 


30 


19 


SBD84, G04a 


G330.0- 


fl5.0 




1200 


63 


180? 


350 


K96 


G332.4- 


-0.4 


2000 


3100 


9 


10 


28 


CDB97, MA86, G04a 


G337.2- 


-0.7 


3250 


15000 


26 


6 


2 


RHS01 


G337.8 


-0.1 




12300 


27 


9x6 


18 


K98b 


G346.6- 


-0.2 




8200 


19 


8 


8 


K98b, D93 


G349.7- 


f0.2 


14000 


14800 


9 


2.5x2 


20 


RM01, G04a 


G352.7- 


-0.1 


2200 


8500 


17 


8x6 


4 


K98a 



a Notes: Many of the radio SNRs have more than one published value for distance 
and age. For these, we either chose the most recent estimates or used an average of 
the available estimates, or the most commonly adopted value. 

b Notcs: Diameters were calculated using from distances together with the angular 
sizes in Green (2006) catalogue. 
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